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Abstract 

We describe a natural particle physics basis for late-time phase transitions in the Uni- 
verse. Such a transition can seed the formation of large scale structure while leaving a 
minimal imprint upon the microwave background anisotropy. The key ingredient is an 
ultra-light pseudo-Nambu-Goldstone boson with an astronomically large (O(kpc-Mpc)) 
Compton wavelength. We analyze the cosmological signatures of and constraints upon 
a wide class of scenarios which do not involve domain walls. In addition to seeding 
structure, coherent ultra-light bosons may also provide unclustered dark matter in a 
spatially flat universe, — 1. 
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I. Introduction 


In a few remarkable instances, modern particle physics theory has predicted the ex- 
istence of new phenomena on macroscopic distance scales. One well-known example 
is the axion, a hypothetical pseudo-Nambu Goldstone boson (pNGB) associated with 
the Peccei-Quinn symmetry, introduced to solve the strong CP problem. 1 Axions arise 
when a global U(1)pq sy mm etry is spontaneously broken by the vacuum expectation 
value of a complex scalar at the scale / a , ($) = / 0 e* a ^®; at this scale, the axion, the 
angular field a around the infinitely degenerate minimum of the potential, is a mass- 
less Nambu-Goldstone boson. QCD instantons explicitly break the global symmetry 
at the scale ~ 100 MeV, generating the axion mass, m a ~ 0(m„f K /f a ). Since its 
couplings and mass are suppressed by inverse powers of /„, the axion is very light and 
very weakly interacting. Nevertheless, it can play an important role in astrophysics and 
cosmology; indeed, astrophysical and cosmological arguments constrain the global sym- - 
metry breaking scale to lie in a narrow window around /„ ~ 10 10 — 10 12 GeV. Thus, the 
axion mass m a ~ 10 -5 eV(10 12 GeV// o ), and its Compton wavelength is macroscopic, 
A a ~ (/ a /10 12 GeV) cm. Although motivated by the strong CP problem, the axion 
is a particular instance of a more general phenomenon, and it portends an important 
lesson: the physical world may contain many new phenomena in the far infrared which 
are not directly accessible in the laboratory, but which may play an important role in 
the development of the early Universe. 

Axions may be generalized to include familons and majorons, 2 as well as more exotic 
objects 3 . Recently, a class of pNGBs closely related to f ami lons (called ‘schizons’), with 
masses of order m# ~ m /erm»on//> has been analyzed in some detail 4 . If one associates 
ITT' fermion with a hypothetical neutrino mass, m„ ~ 0.01 — 1 eV, and / ~ Mqut — M p j ~ 
10 15 — 10 19 GeV, one arrives at a cosmologically interesting scale for the boson Compton 
wavelength, ~ f/™l ~ kiloparsecs - Megaparsecs. This naturally leads to the idea 
of a “late time phase transition (LTPT),” t.e., a vacuum rearrangement occuring at the 



very low temperature T c ~ m v , which may generate structure on a correspondingly large 
scale. 5 If it occurs after decoupling of the cosmic background radiation at Td ec ~ 0.3 
eV, such a late transition opens the possibility of forming large-scale structure without 
imprinting an excessive angular anisotropy 6T/T on the 3 °K microwave background. It 
was originally proposed that “soft” domain walls, with thickness of order ~ //m 9 , 
would form in the discrete symmetry breaking at T c ~ m v and generate non-linear 
density fluctuations, while (it was hoped) avoiding direct imprinting of the walls upon 
the background radiation. 5 

The idea of a late time phase transition is, we believe, more general than the par- 
ticular realizations that have been suggested. However, it is in need of specific detailed 
models and further theoretical refinement before it can be assessed and tested. Pre- 
liminary numerical analysis has suggested that soft domain walls of ultra-light pNGB’s 
may remain potentially problematic, since it has been argued that at least one large 
wall ultimately extends across our Hubble volume 6 , leading to unacceptably large mi- 
crowave background fluctuations 7 . It is unclear, however, whether these simulations are 
ultimately definitive; they do not necessarily have the resolution to see the formation of 
small structures, such as “vacuum bags,” which Widrow and others have argued may 
provide an attractive mechanism for the formation of accretion centers. 8 Moreover, it is 
conceivable that one might still arrive at a viable soft -domain wall scenario with suit- 
able modifications, e.g., if the walls are subject to strong friction due to the surrounding 
medium, or if they divide regions of differing neutrino mass 9 . In this paper, we will not 
consider domain wall scenarios further. 

Recently, Press, Ryden, and Spergel (PRS) have considered an alternative possi- 
bility of a second-order late time phase transition and its consequences for large-scale 
structure, 10 without invoking domain walls. In their model, non-linear fluctuations in 
an evolving soft boson field directly generate large-scale structure (and voids) on scales 
~ 30h -1 Mpc. In addition, the coherent oscillations of the field provide dark matter 
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with critical density, ft^, = 1. The soft bosons are non- relativistic, but their Compton 
wavelength is so large that they do not cluster on galaxy or cluster scales, in agree- 
ment with the fact that the inferred density of matter clustered on these scales is only 
ft ~ 0.1 — 0.2. In this model, the dark matter in galaxies and clusters is thus purely 
baryonic; this is marginally consistent with limits on the baryon density ftj, from big 
bang nucleosynthesis. If viable, this scenario is attractive, since it brings together a 
variety of cosmological problems and solves them in one model. Yet, an important issue 
is whether such a model is reasonable from the viewpoint of particle physics. 

Our only rational guideline in thinking about ultra-low mass particles is the principle 
of “Naturalness.” This is a well-defined operational principle in theoretical elementary 
particle physics, first stated by ‘t Hooft 11 . In this regard, small mass scales must be 
“protected” by symmetries, such that when the small masses are set to zero they cannot 
be generated in any order of perturbation theory, owing to the restrictive symmetry. 
We will not enter here into a general or complete discussion of this mechanism, and how 
to insure its implementation, since the literature of particle theory is infused with this 
principle (see eg., refs. (4, 12) for a more lengthy discussion of naturalness in the schizon 
models and in thermal physics of soft— bosons). In its strongest form, the principle of 
naturalness requires that small mass scales (or large hierarchies) must appear as a 
consequence of some plausible mechanism (in addition to being protected once they 
appear). While the cosmological implications depend upon dynamics and are generally 
insensitive to whether or not a given model Lagrangian has been fine-tuned, the form 
of any given low-energy effective Lagrangian, and its finite temperature corrections 12 , 
will be strongly influenced by the symmetries of the interactions of the full theory. Let 
us summarize the constraint of naturalness in the present context. 

As mentioned above, there can exist a general class of pNGBs (“schizons” 4 ), with 
masses m, ~ m / er m«on//») where the decay constant /, is naturally associated with 
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the GUT or Planck scale, f, ~ 10 15 — 10 19 GeV. Here, the small schizon mass is pro- 
tected by fermionic chiral symmetries or additional discrete symmetries and is therefore 
technically natural. That is, when certain fermion mass terms are set to zero in the 
Lagrangian, the schizon mass goes to zero; the fermion mass terms will not be generated 
in any order of perturbation theory. Familons or schizons are really not significantly 
less compelling than axions; indeed, if one accepts the existence of axions, then it would 
seem undemocratic of nature not to supply pNGB brethren such as familons or schizons. 

From this perspective, the “soft boson” model of PRS, while cosmologically inter- 
esting, is highly unnatural. The Lagrangian is that of a self-interacting complex scalar 
field, with potential V(<f> t<£) = A (<f>^ (f> — v 2 ) 2 and the scale v ~ 10 17 GeV. The mass 
term for the scalar field, m $ — 2A 1 / 2 u, is fine-tuned to be of order m # ~ 2 x 10“ 28 
eV ~ (30 kpc) -1 , by fixing the self-coupling constant to be A ~ 10 -108 . Moreover, 
the field is assumed to have normal strength interactions with other particles. In any 
quantum field-theoretic version of the model, these interactions would lead to an un- 
controllable quadratic divergence of the mass term (and a logarithmic divergence of the 
self-coupling). Thus, to maintain the small mass term one must fine-tune the theory in 
each order of perturbation theory. Also, even if one is willing to accept the unnatural 
fine-tuning of the Lagrangian, the initial condition imposed on the field, <f>i/v 10~ 3 ° 
is ad hoc and fine-tuned. Nevertheless, given the interesting cosmological consequences 
of the PRS scenario, we feel it is worth pursuing more plausible particle physics models 
which can incorporate its attractive features. 

Thus, in this paper we attempt to “naturalize” and amplify the proposal of a late 
time phase transition. We will largely dispense with domain walls and instead follow 
the route of PRS without, however, invoking unnatural fine-tunings of parameters or 
of initial conditions. In Sec. II, we exploit the aforementioned properties of pNGBs to 
construct a class of well-defined LTPT models which are acceptable from the point of 
view of naturalness. In Sec. Ill, we abstract the general features of these models and 
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discuss the resulting cosmological scenarios in detail. We follow with the conclusion, 
in which we also speculate on other possible consequences of an ultra-light pNGB field 
which bec am e dynamical at recent epochs. In the second paper of this series (hereafter, 
paper II), we study structure formation arising from the dynamics of a light pNGB field 
under a wider variety of initial conditions, and we explore concomitant constraints from 
the microwave background anisotropy. 


II. Models 

Above, we argued that an ultralight boson mass scale naturally arises only in models 
which implement symmetries in an appropriate way. In this section, we review two 
classes of toy models with the desired properties: neutrino schizon models and a hidden 
axion model. 

Consider first the low energy effective Lagrangian which contains a neutrino field u: 

C — + VRi^dpVR + Tn 0 v L v R e' p/ f + h.c. (2.1) 

where V(r,l) 0X6 respectively right- and left-handed projections, V(r,l) = (1 ± 7 5 )t// 2 . 
The term proportional to mo can arise from a Yukawa coupling gi>L^R^ + h.c., where 
the complex scalar field $ has a non-zero vacuum expectation value, ($) = fe /V2, 
and m 0 = gf / \/2. This is a familiar chiral Lagrangian, possessing the continuous chiral 
symmetry: 

— * ► e ,ot UL] ur -*• e~' a u R \ <fr -+ <f> + 2 af (2.2) 

Such a theory has several important and well-known properties: (l) it can be embedded 
in a fully renormalizable theory in which a 17(1 )-invari ant complex field develops a 
vacuum expectation value, ($) = /, and <j> is then the residual Nambu-Goldstone boson; 
(2) C is itself renormalizable for a small cut-off A, up to suppressed counterterms of 
order A//; (3) <f> will be identically massless unless terms are introduced which explicitly 


6 



break the chiral symmetry; (4) (f> satisfies “Adler decoupling,” i.e., we may replace v 
everywhere by u': 


v L = u L e t<t,/2f -, 


v'r = VRC t4>/2f ] 


(2-3) 


and we thus see that <f> disappears in the mass term but couples derivatively to the neu- 
trino as 5 #4 ^x/ , 7 5 7 m i/'//' Therefore, for small <f> momentum q tl) <f> emission or absorption 
amplitudes will tend to zero. As a consequence of this decoupling theorem, <f> will not 
mediate a long range 1/r 2 force (we note, however, that the decoupling theorem can be 
violated when the symmetry is broken by a chirally non-invariant mass term). 

Let us now consider explicitly breaking the symmetry. To the Lagrangian of (2.1) 
we may add a small mass term of unknown origin. Usually this comes from some deeper 
symmetry breaking in the theory which breaks the continuous U(l) down to a discrete 
subgroup Zn. For example, let us break U( 1) — > Zi- This implies that <f> — ► <f> 4- mrf 
remains an invariance. So we now have: 

(2.4) 


C = C + K* 


COS 



The physical mass of <f> is determined by shifting f to a local minimum and then ex- 
panding the cosine to quadratic order in <f > . We obtain: 


m# = 2n 2 1 f (2.5) 

and the quartic interaction term A<^ 4 /4!, where 

X = 16* 4 // 4 . (2.6) 

The limit k — ♦ 0 is the symmetry limit of the theory in which we recover the full 
continuous U(l) chiral symmetry. Therefore, k can be naturally small in the technical 
sense; radiative corrections from the full theory, or even for the effective theory, will only 
multiplicatively renormalize «, since when k = 0 the symmetry prevents these effects 
from generating a nonzero « (hence, small /Cj, ore produces a small « r ^normalized «&are 
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in perturbation theory). Of course, as we stated above, the mass term comes from some 
deeper structure in the theory, and when we take the symmetry limit we know that this 
deeper symmetry breaking structure has also gone to its symmetry limit. However, just 
adding the mass term for <f> has one special property: it breaks the symmetry only in 
the mass term of 0 and the decoupling theorem will still hold. 

What kind of deeper structure might give rise to such a mass term for 0? In the case 
of QCD, the proton and neutron are analogues of the v field and the pion is the analogue 
of 0. The deeper structure that breaks the chiral symmetry there is the presence of light 
quark masses, which are not chirally invariant. This leads to the nonzero pion mass, and 
a small non— chirally invariant contribution to the nucleon mass (the cr term). We can 
make an analogy to this situation in the present case by adding an explicit neutrino mass 
term to the Lagrangian that breaks the chiral symmetry. The low energy Lagrangian 
then becomes: 

C =|$ M 0d*0 + VL^d^UL + VRi^d^VR 

7 + e 

Notice that the C/(l ) symmetry is now broken here to the trivial center by the neutrino 
mass terms as well as the cosine term, and therefore only the residual discrete symmetry 
0 — » <f> + 2mrf remains as an invariance. (As a consequence, the argument of the pNGB 
potential is now 0/ f instead of 20//.) 

Now, if k — ♦ 0 we must also set e — ♦ 0 to recover the true sy mm etry limit. However, 
a nonzero k will always be induced at one loop by the presence of a nonzero e and mo. 4 
With a cutoff A < /, we find the induced term: 

Ci-I-P = <M*lf) • (2.8) 

We can freely view this as the origin of the scale k 2 ~ y / moeA. Thus far, this is a 
neutrino version of the schizon model of Hill and Ross 4 . We note, however, that we 


+ (m 0 u L u R e t,f, / f + zvlvr + h.c .) + k 4 cos 
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have ignored the possible effects of CP-violation which can lead to induced Yukawa 
couplings of 4> to Du, and thus to new long-range forces between neutrinos of quasi- 
gravitational strength 4 ’ 13 ; a more general model should allow for a CP-phase on the e 
term as discussed below. 

In this theory we see that the induced scalar mass will be of order: 

~ m 0 e(A 2 // 2 ) ~ m 0 e. (2.9) 

The mass for (f> can now be technically naturally small since we can tune the symme- 
try breaking parameter e to be arbitrarily tiny for large mo; e.g., the observed neutrino 
mass is m„ ~ m 0 ~ 1 eV while ~ (100 Mpc) -1 if e ~ 10~ 8 ° eV. The symmetry 
guarantees that radiative corrections will not change this result. 

This theory will generally have a late-time phase transition at a temperature of order 
the neutrino mass, T ~ mo ~ 1 eV. 12 However, having to input the small parameter e 
is still riot completely satisfactory. Preferable is a scheme in which the scale of ~ 100 
Megaparsec is generated by a strongly natural mechanism that involves only the putative 
scales of particle physics, such as m„ ~ 1 eV and Mqut o r Mpi anc k , in concert with 
some additional symmetry principle. 

In the Lagrangian of eqn.(2.7), we observed the appearance of a (“large”) quadrati- 
cally divergent contribution to the induced mass of (j> , eq.(2.8). Can we somehow reduce 
the degree of divergence of this induced term? Yes, in fact residual symmetries can 
readily control this. Consider first the following schizon Lagrangian invariant under a 
residual Zi discrete symmetry: 

2 ( 2 . 10 ) 

+ (m 0 + ee'*/ f+l 3 )u L UR + (m 0 - ee t4>/f+0 )XLXR + 

where the allowed CP-violating phase, (3, is arbitrary. The continuous 17(1) chiral 
symmetry is broken down to a residual Z 2 discrete symmetry: 

V ~+X\ X^ v \ 0-+0 + 7T/ (2.11) 

9 



If one now computes the induced <f> mass term, one observes that the A 2 term of eqn.(2.8) 
cancels. The leading contribution is now only log-divergent: 

Cx-ioop ~ ^5- ln(A 2 /m§) cos(2 <f>/f + 2(3 ) . (2.12) 

The <(> mas s scale that is now induced is of order 

ttij, ~ m 0 e/f. (2.13) 


In this case, if e ~ mo ~ m„, we retrieve the desired result that ~ m 2 //, so that an 
ultralight schizon emerges without inputting a tiny mass scale. 

These models can be easily generalized to further soften the contribution of fermion 
loops and eliminate the cut-off dependence in the induced <f> mass altogether. Consider 
the £jv-mvariant chiral theory of N neutrinos, 

N - 1 

C '= + (mo + VjL^jR + h.c. (2.14) 

(Hereafter, we suppress possible CP-violating phases (3.) The 17(1) chiral symmetry is 
broken to a residual Zn discrete symmetry: 

uj -» v j+ i ; * vq ; 4> -*■ <j> + Injf/N . (2.15) 


The 1-loop correction is now 



(2.16) 


where 

M 2 = m 2 + e 2 + 2m 0 ecos + -jjr'j • (2.17) 

This respects the discrete symmetry. For N > 2, the sum EjMj is independent of <^; 
thus, the (^-dependent term is independent of the cutoff A, and for N > 2 we can write 


vw) = - E 

j 


Mj 
I67 r 2 


In M) . 


(2.18) 
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In this case, the ^-potential is explicitly calculable. 

In addition to neutrino schizons, one can readily envisage other candidates for “natu- 
ral” low mass elementary particles. For example, consider 14 a hidden, unbroken SU( 2) 
gauge theory (“quasicolor”) containing n'^ massless fermion flavors (“quasi-quarks”), 
which unifies with SU( 3) c at a scale Mgut • By the usual renormalization group, this 
new force will become confining at an energy scale 

/A ~ n \ (33-2n,)/(22-2n',) 

A su ( 2) ^ M GUT ■ ( 2 * 19 ) 

Taking Mgut ~ 10 17 GeV, n/ = 6, n'j = 4, and using A qcD — 0.1 GeV, we obtain 
A5[/(2) « 0.1 eV. If the “quasi-quarks” are massless, then this theory will contain mass- 
less quasi-pions, some diquark “quasi-baryons” with masses of order 2A su( 2 ), and an 
analogue of the rj' with a mass of order A su( 2 ) (due to macroscopic instanton effects). 
These are all “quasi-hadron” boundstates on the scale of Ast/( 2 ), and all of the phe- 
nomena will be natural in the strong sense. A quasi-hadronic phase transition would 
occur when the Universe has a temperature T ~ A st/( 2 )> **e., at a redshift z ~ 400. If 
such a phase transition occurred within the hidden sector of a unified theory, e.g., if 
the SU{ 2) is the low energy remnant of the hidden Ea of superstring models, we would 
only know of its existence through the gravitational effects that it produces. If the 
9 parameter of the 517(2) theory becomes associated with a dynamical field through 
a Peccei-Quinn-like mechanism, i.e., via an anomalous global symmetry with sponta- 
neous symmetry breaking scale f qa , the resulting pNGB, the ‘quaxion’, will obtain a 
mass m qa ~ ^su{ 2 )l /$«• ^ or /90 ~ Mgut ~ 10 17 GeV, we again find m qa ~ 10 -28 eV, 
with a Compton wavelength of order 30 kpc. 

It is also possible that such cosine potentials can be generated nonperturbatively in 
a theory if the associated symmetry has an anomalous current; this happens in QCD 
where instantons nonperturbatively generate a large 77' mass, or a nonvanishing axion 
mass, through the axial U( 1) anomaly. In the present context the prefactor of the 
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cosine potential might be the large scale A 4 ~ / 4 multiplied by an extra “tunneling 
suppression” factor of order exp(— 87r/a), yielding a very small mass for the pNGB. 
The application of ideas like this to soft-boson models has been suggested by Ovrut 
and Thomas 15 ; this may ultimately be a very appealing way of generating soft-boson 
mass scales. 

To study phase transitions in these models, we need to know their behavior at fi- 
nite temperature T. 12 For the neutrino schizon models, this involves integrating over 
a thermal Fermi-Dirac distribution of neutrinos to get the one-loop finite temperature 
effective potential Vt(4>). There are two subtleties here: first, one is interested in study- 
ing effects at T ~ m„ ~ eV, yet the light neutrinos freeze out of thermal equilibrium 
at Tp ~ 1 MeV. Although the neutrinos are no longer in thermal equilibrium, their 
distribution function retains the form of a redshifted thermal distribution so long as 
they are semi-relativistic. (This distribution is characterized by an effective temper- 
ature T e ff — Tpa{tp) I a(t), where a(t) is the cosmic scale factor; one thus computes 
Second, the reactions that keep <f> in thermal equilibrium have rates of order 
T ~ T 3 // 3 . The field will be in thermal equilibrium if this rate is larger than the Hubble 
expansion rate H, where, for a radiation-dominated universe, H ~ T 2 /Mpi anc k- Thus, 
the condition that cf> be in equilibrium, T ^ H, corresponds to T ^ f 2 / Mpj anc fc. Hence, 
for / ~ MguTi <t> decouples thermally very early. However, we emphasize that this 
does not invalidate use of the effective potential. In the neutrino— schizon models, one 
computes the temperature corrections by calculating the value of operators to which <f> 
couples, such as uv, in the appropriate density matrix for the neutrinos. We emphasize 
that this is a coherent or classical field treatment of <f> and uu, corresponding to a clas- 
sical limit of quantum mechanics, and it is not invalidated by the small reaction rates 
for (f> particles to scatter incoherently off of a given v excitation 16 . 
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The relevant temperature corrections have been calculated for the Z at neutrino 
schizon models in ref. 12. In the high temperature limit, T mo, the result is 


A V T (<t>) = 


E 




( 2 . 20 ) 


Comparison with eqn.(2.18) shows that, for the Zi model, the coefficient of the cosine 
potential changes sign at a critical temperature T ~ (mo ± e): the schizon undergoes a 
second order phase transition at a temperature comparable to the neutrino mass. If the 
zero-temperature minimum of the potential is at <f>/ f = 7r, with maxima at <f>/ f = 0,27r, 
then at temperatures T m„ , <{>/ f = ir becomes a local maximum of the potential, 

and the zero-temperature maxima become minima. This is illustrated in Fig. 1. On 
the other hand, for the N > 2 models, because of the ^-independence of the expression 
in the sum of the vacuum contribution of eq.(2.18) and the finite temperature 
contributions of eq.(2.20) the <j> dependence cancels at T ^ mo- In this case, the <f> 
potential becomes asymptotically flat as the temperature is raised, as shown in Fig.2. 
As a result, the transition at T ~ m v in this case is analogous to the transition for 
an axion potential arising from QCD instantons at T ~ A qcd in which the axion 
mass turns on. Note that this is also the expected finite temperature behavior for the 
hidden 517(2) quaxion model outlined above. There is, however, a qualitative difference 
between the Zw >2 and quaxion models: for the schizon models, the potential turns 
off only logarithmically as the temperature is raised, while for quaxions the instanton- 
induced potential is suppressed roughly as an inverse power of the temperature. For 
axion models, it is well known that the temperature dependence of the mass plays an 
important role in determining the cosmological bound on / due to the density of coherent 

i 

axion oscillations. For the schizon models, however, this temperature dependence is 
small: when the schizon starts oscillating, its mass differs only logarithmically from its 
zero temperature value. The cosmological implications of these different behaviors will 
be discussed in the next section. 
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In the model of Ovrut and Thomas 15 , it appears that there is no temperature 
dependence of the potential for T & /; in this case, the “phase transition” will simply 
involve the ‘unfreezing’ of the boson field when the Hubble expansion parameter becomes 
of order the <^-mass. If, however, the tunneling effects which generate the potential are 
associated with a mass scale M , then we would generally expect the potential to turn 
on as a power of temperature at T ~ M , as in the (qu)axion model. 


III. Cosmology and Late Time Phase Transitions 

Having discussed a variety of underlying particle physics models for late phase transi- 
tions, we now turn to their cosmological implications. We will examine consequences 
that are relatively insensitive to the details of particular models. Thus, we assume the 
existence of a generic pNGB with a phenomenological Lagrangian given as a function 
of temperature: 

LpNOB = \dp4>a*<t> + VtW , (3.1a) 

where the finite-temperature potential is 

V T {<j>) = + c(T)[cos(*//) - 1]) . (3.16) 

Here, (f) is a real pNGB scalar field, / is the scale of global spontaneous symmetry 
breaking, assumed to be of order the grand unification or Planck scale, / ~ 10 15 — 10 19 
GeV, and n is an explicit symmetry breaking scale associated, e.g., with a light neutrino 
mass or a hidden strong interaction scale, taken to be of order 10 -2 — 1.0 eV. By suitable 
rescaling of fx we can set c(T = 0) = 1, and we have inserted a constant term in the 
potential to insure a vanishing cosmological constant in the usual way. 

From the discussion of particle physics models in Section II, we extract three broad 
classes of models for the temperature-dependence of the potential: 
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I. Models in which c(T) ~ 0 at high temperatures, T T c ~ /z, and the potential 
turns on approximately logarithmically at a critical temperature T ~ T c ~ /z; for 
these models, c{T) ~ 1 for T T c . The Zn schizon models with N > 2 are in this 
category. 

II. Models in which c(T) is a slowly varying function of temperature that changes sign 
at a critical temperature T c ~ /z, defined by c(T c ) = 0; for example, in the Z 2 
schizon model, c(T) ~ ln(T c /T), where T c ~ fj. ~ m„. For models of this type, the 
Universe goes through a second-order phase transition at the critical temperature 
T e . For T ^ T c , we will generally assume |c(T)| ~ 1. 

III. Models in which c(T) ~ 0 for high temperatures, T T c , and the potential turns on 
roughly as a power law in temperature near a critical temperature, c[T ) ~ (T c /T) n , 
n ~ a few (for T ^ T c ). For these models, c(T) ~ 1 for T ^ T c . For example, 
the quaxion model, where the instant on-induced potential is suppressed at high 
temperature, is of this type. 

We will not explicitly consider a fourth logical possibility, that c(T) is constant 
for all T, as in the Ovrut-Thomas model with a high energy tunneling scale. This 
possibility is implicitly included in some of the cases discussed below. 

Abstracting from the models of the preceding section, we generally expect the critical 
temperature to be of order the explicit breaking scale, T c ~ /z; we will define 

f = T c /n = 0(1) (3.4) 

as a third parameter of the model, which is naturally of order unity. We note that, in 
some cases, the potential term in eqn.(3.1) may actually be absent until the Universe 
cools to the electroweak symmetry breaking scale, T ~ 100 GeV. Indeed, in the schizon 
models this is necessarily so since the potential term arises from quark or lepton masses, 
which do not appear until electroweak breaking. 
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With the generic temperature-dependence of the potential in hand, we now discuss 
the cosmological evolution of the scalar field in broad outline. We focus here upon the 
spatially homogeneous, zero-momentum mode of the field, 4>{t) = (<f>(x, t)), where the 
brackets here denote spatial averaging. In paper II, we study the evolution of spatial 
fluctuations in the field; here, we are implicitly assuming that these fluctuations do not 
strongly perturb the zero mode. This is certainly the case if the fluctuation amplitude is 
small compared to 0 as would be expected, e.g., after inflation if the reheat temperature 
Trh < f- in this case, aside from inflation-induced quantum fluctuations, the field will 
be homogeneous over many present Hubble volumes. On the other hand, if inflation 
did not take place, or if Trh > f so that the global symmetry is broken again after 
inflation, we generically expect large spatial gradients (“Kibble gradients”) in the field 
due to the fact that <f>/ f is uncorrelated on scales larger than the Hubble radius when 
the transition at T ~ / occurs. In this case, it may not make sense to talk of a zero- 
momentum mode? 7 , but , at least in the absence of light boson production by topological 
defects, the zero-mode treatment does lead to an estimate of the scalar energy density 
which should be accurate in order of magnitude, since the (neglected) gradient energy 
term is comparable to the potential energy term 18 . 

We assume that at temperatures / <; T >> T c , ^ is a classical field expectation 
value, randomly placed on its potential. The scalar equation of motion is 

% + 3H<f> + —~rr^ — 0 > ( 3 . 5 ) 

d<p 

where the Hubble parameter is given by H 2 = {Bir/Z m^p for a spatially flat universe. 
At these high temperatures, the Compton wavelength of the field is much larger than the 
Hubble radius, so the potential term is negligible compared to the Hubble damping term 
(4> is effectively massless). In this limit, the solution is time-independent, 4>{i) = 4 >\ , 
i.e., the field is frozen to its initial value. 
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As the temperature drops below T r , defined by the point where 


|mJ(T r )l = \#V T JiP a 9S 2 (T r ) , (3.6) 

the <j> field becomes free to roll down the potential, modulo Hubble damping. There 
will then result spatially coherent oscillations about the potential minimum; once the 
amplitude is sufficiently small that the oscillations axe approximately harmonic, the 
0 stress tensor is that of a nonrelativistic particle (pressureless dust). The coherent cf> 
oscillations may currently dominate the mass density of the universe, providing (unclus- 
tered) dark matter, but they must satisfy f 2^ < 1. This condition provides a constraint 
on the parameter space of /, /x, and £ (see below). 

We now discuss the cosmic evolution of the scalar field for the three classes of models 
in more depth. 


IIL1 MODEL I 


For these models, the potential is flat at high temperatures and turns on rapidly at 
the critical temperature T c (see Fig. 2). We shall approximate their behavior by letting 

c(T) = 0 , T>T C 
= 1 , T <T C 


(3.7) 


The model is described by 3 parameters, /, /x, and £ (or T c ), and it is useful to define 3 
important characteristic temperatures: T c is the critical temperature as defined above 
(i.e., when the potential goes from being flat to curved); T r is the temperature, defined 
above, when the Hubble damping ceases and the field begins to oscillate in the potential; 
and Td is the temperature when the energy density of the scalar field begins to dominate 
over baryons (if ever). By definition, T r < T c . 

We will assume that the initial value of the field, <^>i, is not extremely close to 
the maximum or minimum of the zero-temperature potential. For a typical initial 
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value of the field, the potential energy density at the critical temperature is of order 
Vx e {4> i) ~ /i 4 . It is convenient to define a parameter a which characterizes the ratio of 
the scalar energy density to the baryon density, 


a (T) = 



(3.8) 


Thus, the temperature Td is defined by a(T<f) — 1. Since ps(T c ) = ^lBPerit(T c /To) 3 , 
and using T c = £/i, we have 


a 


(Tc) = 


mT 0 3 


flBPcritt 3 


= 8.5£ 


-3 


/ M ( 0-02 
VeW \Sl B h? 


(3.9) 


where, in the last equality, we have used the present microwave background temperature, 
3o = 2.4 x 10 -4 eV, and the critical density p cr it — 8.1 X 10 -11 /i 2 eV 4 . We also define 
the ratio 

(SI), • (3 - 10) 

Thus, the temperature T r , at which the field starts oscillating, is defined by /3{T r ) = 1. 

It is convenient to divide the discussion of this model into two sub-cases: 1) a(T c ) < 
1, i.e., p. < (^ 3 /8.5)(nsh 2 /0.02) eV: baryons dominate at T — T c and T c > Td\ 2) 
a(T c ) >1:0 dominates at T c , and Td > T c . We consider these two cases separately. 


III.1.1 Case 1: a(T c ) < 1 

Since baryons dominate at T e , we must use the baryon density to determine the 
expansion rate in the denominator of /3 (eqn.3.10); for a typical field value on the 
potential, we have |m|| ~ /* 4 // 2 , and so 

«r.) - «(r.) (j=jt) . (3 - u ) 

For this case, since a(T c ) < 1, we have (3(T C ) < m 2 pi /2A-Kf 2 . Again, we subdivide the 
possibilities: 
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la) / > mpij s/TAk = 1.4 x 10 18 GeV: in this case, we are guaranteed to have 
/?(T C ) < 1, so that T r < T c : once the potential ‘turns on’, the field remains frozen 
to its initial value for a while, because of Hubble damping. As a result, the scalar 
energy density remains approximately constant (instead of redshifting), acting as a 
cosmological constant during this phase. Since the baryon density redshifts, the ratio 
a(T) grows as a ~ T -3 , giving us the more general relation /?(T) = a(T)m 2 pi /24irf 2 , so 
that 0(T) < a(T). Thus a(T) reaches unity before /3(T) does: <£ is frozen to its initial 
value until after it dominates the energy density, i.e., we have the relation T c > Td > T r . 
This can be seen from evaluating 



(3.12) 


while the critical temperature is given by 


Thus Td/T c = a 1 / 3 (T c ) < 1. The evolution of the scalar and baryon densities is shown 
schematically in Fig. 3. 

Once the temperature drops below Td, the scalar energy density dominates the 
expansion rate, and the expression for (3{T ) changes. Initially, at T just below Td, the 
field is still frozen at its initial value, so the density V(<£) ~ n* ~ constant. In this 
case, we find (3(T & Td) = Since / > mpi/y/2Ar, this implies (3{T ^ Td) < 1, 

so the field remains approximately frozen. In actuality, however, the field does evolve 
toward the minimum, but on a timescale slower than the expansion time H ~ 1 . When 
/3(T) ^ 1, the field slowly rolls down the potential with a characteristic growth rate 
^ ~ exp(m^t/3H') ~ exp(3 Ht(3). One can think of this as a brief period of slow rollover 
quasi-inflation beginning a t T = Td and ending when the field begins to roll rapidly 


T c ( M 4 

To \Cl(Tc)ftBPcrit 


1/3 


(3.13) 
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down the potential at T = T r . As we shall see, however, this epoch must be very brief, 
so that the universe never enters a full-blown exponential de Sitter inflation. 

Witfi a simple approximation, we can analyze the slow rollover epoch: we assume 
that the baryon density is negligible compared to the scalar density at T <T d and that 
the scalar energy density remains constant from T = T d down to T = T r , below which 
it subsequently scales as non-relativistic matter. As a result, we have, at the present 
epoch, 

^ = (**) cz(T r /T d ) 3 = e~ 3N - (3.14) 

where N e = ln(T<i/T r ) is the number of e-folds of the scale factor during the quasi- 
inflationary epoch between T d and T r (note that we are not assuming exponential ex- 
pansion during this phase). Since big bang nucleosynthesis suggests f lg ~ 0.01 — 0.21, 
we require 

lV e .~(l±0.5) + ^ln(ft*) (3.15) 

e-folds of quasi-inflation, i.e., the scalar field should slow-roll only for about 1 expansion 
time. 

The slow roll-over of a scalar field in a potential of the form (3.1) has been analyzed 
recently in the context of a model of inflation 19 , and we can apply those results here. 
The slow rollover phase ends, and the field begins oscillating, when the field reaches a 
value given implicitly by \V'($ r )mpi/V(<fi r ) j = \/487r, or 


sin (0 r //) V 487T 


(3.16) 


1-t-cos {4>r/f) rnpi 

Note that this generally happens while jV"(0)( ^ 9 AT 2 , i.e., slightly before the field 
begins oscillating according to our criterion (/3(T) = 1) above. For example, for / = 
mpj, we have <£ r // = 2.98, while for / = mpifyj 247T, <j> r / f = 1.9. If the field begins (at 
T a T d ) at a value 4> = <^j, the number of slow roll e- folds before it begins to oscillate is 


- - I67 r/ 2 

N e ($i ,<£r,/) = s— In 


m 


Pi 


sin(<ft r /2/) 
sin(<£i/2/) J 


(3.17) 
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To achieve N e = 0.5 — 1.5 requires, for example, f>\j f = 2.63 — 2.82 if / = mp\ and 
4>i/ f — 0.17 — 0.79 if / = mpij \/ 2 Note that these values are only indicative: for 
N e ~ 1, the approximation we have adopted that p(<f>) pp is not well justified. 

lb) / < mpi/ \/2An\ In this case, from eqn.(3.11), we have (3 (T) > ct(T), so that 0 
starts to oscillate before it dominates the energy density. Once it begins oscillating, p<f, 
scales like non-relativistic matter, i.e., like pp, so it never dominates the energy density, 
P<t> < PB for all T. (See Fig. 4). In this case, (f> cannot be the dark matter, but it can 
still play an important role in structure formation (see paper II). 

For this case, the parameter space again naturally splits into two regimes: lbi) either 
the oscillations occur after the phase transition, T P < T e , which happens if 0(T C ) < 1, 
or (lbii) the field starts oscillating as soon as the potential turns on, T r = T c , which 
corresponds to (3{T C ) > 1. From eqns. (3.9) and (3.11), in the first regime (lbi), 


24 t ^ T 
8.5 rn 2 pi 0.02 ’ r c 


and the oscillations begin at the temperature 

T r = To ( 


p. A m 2 P i \ 


\2Air f 2 ^Ippcrit 


7 .) 


(3.18) 


(3.19) 


The second regime (lbii) corresponds to 




m 


(3.20) 


Pi 


In this case, the ratio of the Hubble radius to the <f> Compton wavelength when it begins 
oscillating is approximately N r ~ /? 1</2 (T C ), or 


3 mpi / _M 1/a / 0-02 \ a/a 
y/24rc£ 3 f 'eV/ VOsh 2 ; 


(3.21) 


which may be substantially larger than one. This is in contrast to the usual case, e.g., for 
the axion, which starts oscillating when N r = m^/ZH = H~ 1 /3A^ ~ 1. This difference 
has important ramifications when we consider structure formation in these models. 
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Briefly, in the scenario of ref.(10), the scalar field generates non-linear perturbations 
on the scale of its Compton wavelength; unless N r > 1 or the transition occurs before 
recombination, this may lead to unacceptably large microwave background anisotropies 
(see paper II). 

Ili.1.2 ’Case 2: a(T c ) > 1 

In this case, the scalar field dominates the energy density at T c , and Td > T c > T r . 
From the discussion below eqn.(3.13), before the field starts redshifting we have 0(T) 
rrip l /24irf 2 . We therefore split the discussion again into two subcases. 

2a) / > rn.pi/ ■%/247t : in this case, /3(T C ) < 1, so the field is frozen at T c ( t.e ., 
T c > T r ) and may undergo a brief quasi-inflationary slow-rollover. This is similar to 
case (la). Inflation actually begins at T ~ Td given by eqn.(3.12). However, in the range 
Td> T > T c , the potential is flat, and the field does not roll classically. The number of 
e-folds of inflation between T c and T r is constrained by the requirement, analogous to 
eqns.(3.14) and (3.15), that 

MT r ) = a(T e )e JN -< T - r -> = (3.22) 

where a(T c ) is given by eqn.(3.9). Since, for case 2, a(T c ) > 1, this yields the require- 
ment that N e (T c ,T r ) < (l/Sjl^fl^ 1 ) < 0.5 — 1.5; again, the number of inflationary 
e-folds must be small. For fixed \i and £, eqn.(3.22) fixes the number of e-folds re- 
quired to have the proper ratio of scalar to baryon density today; from eqns.(3.16) and 
(3.17), this fixes the required initial value of the scalar field (j>\. Since the number of 
e-folds is small, the initial value of the field must in general be reasonably large, i.e., 
not fine-timed very close to the origin. 

2b) / < mpi/\ / 24ir : here, (3{T C ) > 1, so the field starts oscillating at the critical 
temperature, T r = T c . In this case, the ratio of the scalar to baryon densities is pre- 
served, a(To) = a(T’c)* For < 1, we thus require a(T c ) < 1 /ftp; using eqn.(3.9), this 
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yields the constraint 

A* < = l-^so^eV , n 0 < 1 (3.23) 

- i o 

where the present value of the Hubble parameter is if 0 = 50 /i 5 o km/sec/Mpc, and 
observations indicate 1 < hso < 2. (We use the parameter hso because its default 
value, h SQ = 1, yields marginal agreement with globular cluster ages if 0 = 1.) For 
fixed physical scale /*, we may think of this as determining the required value of (. For 
H = CJ(eV), the scalar field provides closure density if £ = C?(l). Since we expect £ to 
be of order unity, this shows that the natural explicit symmetry breaking scale for a 
coherently oscillating field that dominates the present energy density is in the eV range. 
That this is far below the QCD scale, A qcd ~ 100 MeV, that sets the scale for axion 
oscillations, does not contradict the fact that axions may also dominate the density (if 
/ ~ 10 13 GeV), because the temperature dependence of the axion potential is different 
(see model III below). 

For closure density bosons (ti# = 1), using £(/i) from eqn.(3.23) we find the critical 
temperature 

T c = (400 = 0.9ft s - !/I eV (^) 4/1 , (fl+ = l) (3.24) 

and the redshift of the transition, 

1 + = 3.7 4/3 , (n* = l) (3.25) 

If we require the transition to occur after recombination, i.e., at z c < 1000, and before 
the first known quasars, i.e., z c > 5, eqn.(3.25) gives 

8 x 10“ 3 eV < [ihy j 1/2 < 0.4 eV (3.26) 

This also corresponds to the range 

0.2 < (hlf < 0.7 (3.27) 
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(For completeness, for case 2a above, the corresponding constraint is 5 < z r < 10 3 , 
which yields a range for y identical to eqn.(3.26).) For reference, we note that the ratio 
of the Hubble radius to the scalar Compton wavelength at the onset of oscillations is 
given in this case by 

K a 0 1/J( T C ) = (3.28) 

For case 2b, this is larger than one. 

The parameter space for the different cases of Model I is displayed in Fig.5. 


III. 2 MODEL II 


For these models, the high temperature potential does not vanish, but is inverted 
from the low temperature potential, as in Fig. 1. We approximate the behavior of this 
class of models by choosing the temperature-dependent coefficient in eqn.(3.1) to have 


the form 


c(T) = -1 , T > T c 

= 0 , T = T C (3.29) 

= i , r<r c 

The analysis is similar to that of Model I, the primary difference being that in some cases 
the field may now oscillate before the critical temperature is reached. We must now 
define two temperatures: Tp ^ (Tp 2 ^) is the temperature when the field starts oscillating 
in the high-temperature (low-temperature) potential. Clearly, Tp 1 ^ >T C > Tj?\ 

We divide the analysis into the same subcases as for Model I. Then, for cases la, lbi, 
and 2a, the behavior of Model II is identical to Model I, because the field is always frozen 
in the high temperature potential (/ 3(T > T c ) < 1); for these cases, the constraints are 
summarized in Fig. 5. We now turn to the cases where the field does oscillate in the 
high temperature potential. 
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lbii) For p in the range given by eqn.(3.20), the field begins oscillating in the 
high temperature potential, corresponding to > T c = Tr 2 \ The onset of these 

oscillations occurs when (3(Tr = 1 , which implies 

/i 4 mpj 

fWcr»t247T/ 2 

As in Model I, for this case the scalar oscillations never dominate the energy density of 
the universe. Nevertheless, it is of interest to study the behavior of the field. 

When the field starts oscillating at T^\ typically there will be rough equipartition 
between the kinetic and potential energy of the field, and the energy density decreases 
with the expansion. After a short time, the kinetic energy has redshifted to become 
negligible compared to the potential energy, and the energy momentum tensor of the 
field approximates that of a vacuum state, with vacuum energy p vac ~ 2/x 4 . This does 
not lead to inflation, however, because by construction the scalar field never dominates 
the total energy density for this case. Defining p $ — p# — p V ac, it is easy to show that 
p<f> redshifts like non- relativistic matter, p<p ~ i2 -3 . This implies that the oscillation 
amplitude decays as 

- ! T \ 3/2 

m = (t?) (3 ' 31) 

during the high temperature phase. Using eqns.(3.13) and (3.30), the amplitude at the 

critical temperature is then 



) 


1/3 


(3.30) 


}(TP =r „2 (Tc) = ^n in /evn^ y /! 


/eV n B h 2 
v p 0.02 


<pi ’’ v_fc/ mpi 3 ( 3 - 32 ^ 

In this case, zz T c , and the ratio of the Hubble radius to the Compton wavelength 
at the onset of the low temperature oscillations is just the factor by which the amplitude 
has been damped in the high temperature phase, 


N r = = (^1) 


(3.33) 
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2b) As for Model I, this is the most interesting case, since the scalar field can dom- 
inate the energy density. As above, we have > Td > T c — Tr 2 ^; the situation 
is qualitatively similar to Model I (2b), except that in this case due to the high tem- 
perature damped oscillations, the field may be localized near the origin at the critical 
temperature, 0(T C ) /. This factor is straightforwardly estimated: 

(3,4) 

01 \ m Pl a { T c)/ ™Pl \“<*>/ 

where in the second equality we have used a(T c ) = In this case, the ratio of 


the Hubble radius to the Compton wavelength at T c = Tr is 


N r = = 


mpi 


01 ( 


1/2 


(3.35) 


\/24ir/ <j>{T c ) 

For some applications in paper II, we will be interested in the case where the high 
temperature (in addition to the low temperature) oscillations begin after recombination, 
i.e., zi ^ == TW/Tq < lO 3 .; Using eqn.(3.30), this happens if 

1/2 


/x < 0.2 eV 


fV2&h 50 \ f n B \ 17 4 
V mp, ) VO-0 6j 


(3.36) 


III.3 MODEL III: QUAXIONS 


For this class of models, we take the coefficient in eqn.(3.1) to be 

<T) = (^) , T > T< 


(3.37) 


= 1 , T < T c 


From the effective Lagrangian viewpoint, the parameter n > 0 corresponds to a fourth 
parameter of the model (compared to three for models I and II); however, it is deter- 
mined by microphysics. In this case, the potential is flat at T T c and gradually turns 
on, reaching full strength (and its zero-temperature value) at T c . For field values away 
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from the maxima of the potential, it is convenient to expand Vr(<f>) around the potential 
minimum, = 4> — tt/, 

1 (3.38a) 


VtW = f{T) + - 2 Tn\(T)<f 


where 


f(T) = 2(X 4 [l - (§) 
m*(T) = £ (f ) 


n/a 


T>T C 


(3.386) 


We define the parameters a(T) and /3 (T) as before. In this case, like the axion, the 
field can start oscillating at T r > T e , and it is convenient to divide the cases slightly 
differently than above 20 : 1) a(T r ) < 1, i.e., a baryon-dominated universe when the 
oscillations begin, and (2) a(T r ) > 1, in which case the scalar field dominates the 
energy density when it begins to oscillate. For a typical initial field value, p<t>(T r ) ~ /x 4 , 
so that 


“ (r ' )= te),= 


n B Pcrit(Tr/T 0 y 


(3.39) 


III.3.1 Case 1: a(T r ) < 1 

In this case, p<j,(T r ) < p B {T r ). The temperature T r when the field starts oscillating is 
determined by (3(T r ) = m^(T r )/9if 2 (T r ) = 1, where the temperature-dependent mass is 
given by eqn.(3.38b) and the expansion rate at T r is fixed by the baryon energy density. 
Straightforward manipulation gives 

ll/(3+n) 


Tr 

T c 


rripi 8 


.5 / p, \ ( 0.02 \ 

vTv/ VR^V. 


L24tt/ 2 £* 

' m 2 Pl 8.5 f f 0.02 X ] 1/3 
,24: r/ 2 ? \eVJ \n B h 2 ) . 


, T r > T c 
Tr < T c 


(3.40) 


From (3.39) and (3.40), we have 




m 


PI 


(3.41) 
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so the condition a(T r ) < 1 requires 



(3.42) 


As before, we subdivide this case into two subclasses depending on the value of /: 

la) / > mpt/ >/24ir : In this case, from (3.42) we immediately have T r < T c . Since 
the field does not begin oscillating until the potential reaches its zero-temperature form, 
the high-temperature behavior of the potential is irrelevant in this instance. This is 
therefore identical to Model I, case la. 

lb) / < mpi/ y/2Air : in this case, the ratio T r /T c may be smaller or larger than one, 
depending on the value of p. From eqn.(3.40), the condition T r = T c clearly corresponds 
to 


£ 3 eV fCl B h 2 \ (2±Ttf 2 

fj, _ Peru - g 5 ^ Q Q2 J ^ m 2 pi 


T r = T c 


(3.43) 


(Compare eqns.(3.18) and (3.20).) We thus subdivide this case into two subclasses:(lbi) 
p < Peril, which implies T r < T c , and (lbii) p > p C rit, » e., T r > T c . 

lbi) p < Pcrit- this is identical to case (Ibi) of Model I. Again, the high-temperature 
nature of the potential plays no role, since T r < T c . 

lbii) p > Pcrit ‘ in this case, the field starts oscillating while the potential is still 
evolving. The condition a(T r ) < 1 restricts p to the range 


( m Pl \ 
\247r/V 


(n+3)/n 


> 


JL 

Pcrit 


>1 , T r > T c 


(3.44) 


Although the scalar field never dominates the energy density in this case, it is of interest 
to study the damping of its oscillations in the evolving potential. As before, if we 
subtract the instantaneous zero-point energy by defining p# = p# — f(T), then for the 
harmonically oscillating field we find 

p* + 3 Hp* ^1 - = 0 (3.45) 
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which implies 


m 


P* « 


*(?) 

i2 3 


(3.46) 


Since p $ ~ m^(T)^ 2 , the oscillation amplitude decays as 0 ~ m^ l ^ 2 {T)R 3 / 2 . Using 
eqn.(3.38b), this gives the amplitude at T c , 

(n+6)/4 


0(£) 

01 


=(§) 


(3.47) 


and the relative energy density 


«(T e ) _(T C \ 

*(Tr) \T r J 


n/2 


(3.48) 


where the ratio in parenthesis is given by the first eqn.(3.40) and we have assumed 
0i ~ /• As in the axion case, the commencement of 4> oscillations in the evolving 
potential reduces the relative scalar energy density. 


III.3.2 Case 2: a(T r ) > 1 


In this case, the scalar field already dominates the energy density when it begins to 
oscillate, Tj > T r . From eqns.(3.10) and (3.38), for Tj > T > T r , we find 


/3(T) = 


rrl 'pi I Tc 


24t r/ 2 V T ' ’ T > Tc 


m 


(3.49) 


Pi 


24tt/ 2 


T < T c 


where we have used the fact that p<f> ~ /x 4 for T > T r . 

As usual, this motivates us to consider two subcases: 

2a) / > mpif : This case is similar to Model I, case (2a). That is, in order 
to end up with an acceptable value of f the scalar field must undergo a very 
brief period of slow rolling quasi-inflation. The only difference is that, in the present 
case, the field can in principle start evolving classically at T > T c . However, from 
eqn.(3.49), the rollover rate in this high-temperature regime is exponentially suppressed, 
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since ~ exp(3 Hi/3) and (3 ~ T~ n . As a result, in practice, the initial value(s) of 
required to give ^ 1/(0. 01 — 0.2) in this case will be similar to that in Model I. 

2b) / < vnpi/yj 247T : in this case, from eqn.(3.49), we have 

r, ( m p, \ 1/n 

T c V24,r/2 ; 


> 1 


(3.50) 


and the field starts oscillating above the critical temperature. In the temperature range 
Tr> T > T c , the field oscillates in the evolving potential, and its energy density is 
approximately (Cf. eqn.(3.38)) 


H = f(T) + mJ(T 


(3.51) 




(3.52) 


The oscillation amplitude (j> decays according to eqn.(3.47); again assuming a typical 
initial amplitude ~ / and using eqn.(3.50), we find 

p2 \ ( 8 + n )/ 2n 

*PI 

Now the analysis follows along the lines leading to eqn.(3.23): we form the ratio a(T e ) = 
P<I>[T C )I pb{T c ) and use the fact that this ratio is conserved for T < T c , *.e., a(T 0 ) = 
a(T c ). Then setting < 1, i.e., a(T c ) < l/flg, yields the constraint 

< ^4^ f = 1.4fcl„£ J eV f 


(6+n)/2n / 2 \ (S+n)/2n 

.. , <» PcHt / ) “ 

M - T 0 3 V247T/2 y “ ' V24t r/2 / 


(ft* < 1) (3.53) 


Unlike eqn.(3.23), the cosmological constraint on the scale p here depends on /, as in 
the usual axion case; in particular, as / is decreased below the Planck scale, the upper 
bound on p is relaxed 21 , depending on the index n. 

Assuming the scalar field saturates closure density determines the parameter £ and, 
analogously to eqns.(3.24) and (3.25), we have the critical temperature 


Tc - fiM - 0.9 h 50 2/3 eV (^) ( m / ) 


2 \±+i 


, (ft* = l) (3.54) 
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and the redshift of the transition, 


T 

1 + 2e = -f = 3.7 x 10 4 /i 
To 


- 2/3 
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4 /3 /24tt/ 2 
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(3.55) 


IV. Conclusions 

We have studied a variety of particle physics models which can potentially generate the 
observed large-scale structure in the context of late time phase transitions involving 
ultra-low mass bosons. Some of these models can simultaneously explain the dark- 
matter in the Universe. These models are distinguished in containing a fundamental 
large distance scale, d, which emerges as a hierarchical ratio of microscopic scales, e.g., 
d ~ MauT/mi' Symmetry principles, abstracted and borrowed from experience with 
QCD, allow us to build underlying theories of such a phase transition at a very low 
energy scale, T c ^ eV, that are technically and even strongly natural. This low energy 
scale corresponds to a moderately recent cosmological epoch, z c <> 1000, and such 
models may therefore have other potentially striking observational signatures, as in the 
scenarios for large-scale periodic structure discussed in ref.(21). 

The common element of the models considered here is the existence of an ultra-light 
pseudo-Nambu-Goldstone boson, a lighter cousin of the conventional QCD axion, which 
becomes dynamical at this late epoch. We have analyzed the cosmic evolution of this 
field in three broad classes of models and have identified the regions of parameter space 
for which it can make a significant contribution to the energy density of the Universe. 
Unlike the conventional axion, an ultra-light boson associated with global symmetry 
breaking at the GUT scale, / ~ 10 18 GeV, can provide ~ 1 without invoking special 
initial conditions. We note that such a field is generically very weakly coupled and thus 
difficult to detect. 
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For completeness, we close with several comments about other possible implications 
of, and constraints upon, these models. The first have to do with topological defects. If 
inflation does not take place, or if it does occur but with reheat temperature Trh > /> 
the boson field is not expected to be homogeneous over scales larger than the Hubble 
radius at the critical temperature T c . In this case, the symmetry breaking at T c ~ \jl 
can lead to the formation of domain walls, which may endanger the isotropy of the 
microwave background. As in axion models, there are a number of ways to avoid 
this domain wall problem 1 . For example, one can have additional terms in the scalar 
potential which explicitly break the residual discrete symmetry 22 (bias the potential), 
lifting the vacuum degeneracy; this effectively destroys the topological stability of the 
walls, driving them into the false vacuum regions. As long as the vacuum asymmetry is 
large enough, the walls can disappear before they do any damage. 21 Alternatively, if the 
broken global symmetry is CT'(l), then cosmic strings form in the transition at T ~ 
these become the boundaries of the walls which form later at T ~ /x. In models with 
N = 1 minimum of the boson potential, there is one wall per string, and the wall-string 
system destroys itself before it dominates the energy density. 23 We note that in this 
case the strings radiate bosons until the boson mass turns on at T P , and the resulting 
cosmic energy density of string-radiated bosons may be larger than that calculated 
from ‘misalignment’ production in section III by a factor ~ ln(/< r ) ~ 100. 17 In order 
for the gravitational effects of the strings not to perturb the microwave background, the 
symmetry breaking scale is then restricted to roughly / ^ 10 1T GeV. 

A final speculative possibility concerns the interesting cosmological role of an ultra- 
light field which has not yet become dynamical (is still frozen) or only became dynamical 
at moderate redshift (say, z ~ 8 — 10). The origin of the smallness of the cosmological 
constant is still shrouded in mystery, but let us suppose that there is some mechanism 
which sets the ultimate true vacuum energy density of the universe to zero. Then, at 
late times, classically we expect the vacuum energy to be dominated by the lightest field 
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which has not yet evolved to its ground state. If there is an ultralight pseudo- Goldstone 
boson which has yet to start oscillating, i.e., if it satisfies T r To = 2.4 x 10 -4 eV, 
then it will act as a (temporary) cosmological constant, with present vacuum energy 
of order p vae ~ /a 4 . If p. ~ 3 x 10 -3 eV, then O cac ~ 1. This may have potential 
benefits such as increasing the age of the universe (for fixed value of the Hubble pa- 
rameter) and boosting the large-scale transfer function for density fluctuations. The 
simplest way to implement this is to have T c < To, which requires £ = T c //z < 0.08; one 
could relax this bound somewhat if the global symmetry breaking scale f % mpi/ y/2Arc. 
Alternatively, if the field became dynamical at moderate redshift, it could lead to an 
epoch of cosmological loitering, in which the scale factor of the universe passes through 
a moderate Eddington-Lemaitre-like coasting phase, allowing extra time for the growth 
of large-scale structure 24 . 
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FIGURE CAPTIONS 


Figure 1: The temperature-dependent potential Vr((f)) = V{<p ) + A Vt(4>) for the 
pseudo-Nambu-Goldstone field in the Z 2 schizon model. The field exhibits a second 
order phase transition at a temperature T ~ mo; this is characteristic of models we 
classify as Type II (from ref.(12)). 

Figure 2 : The temperature-dependent potential for the scalar field in a Z 3 schizon 
model; this is typical of Zn models with N > 2: the potential is flat at high temperature, 
developing a minimum at T ~ mo- We classify these as Type I models (from ref. (12)). 

Figure 3: The evolution of the scalar and baryon energy densities p<f, and pg as a 
function of the scale, factor a(t) for model I, case la. In this case, the scalar field 
presently dominates the energy density. 

Figure 4: The evolution of p$ and ps for model I, case lb. Here, the scalar field never 
dominates the energy density of the universe. 

Figure 5: The parameter space for models I and II, showing the different cases in the 
plane / (GeV) vs. p/( 3 (eV). 
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